Increasing stability in an inverse problem for 
the acoustic equation 

Sei Nagayasu* Gunther Uhlmann^ Jenn-Nan Wang^ 



Abstract 

In this work we study the inverse boundary value problem of determin- 
ing the refractive index in the acoustic equation. It is known that this inverse 
problem is ill-posed. Nonetheless, we show that the ill-posedness decreases 
when we increase the frequency and the stability estimate changes from log- 
arithmic type for low frequencies to a Lipschitz estimate for large frequen- 
cies. 

1 Introduction 

In this paper we study the issue of stability for determining the refractive index in 
the acoustic equation by boundary measurements. It is well known that this in- 
verse problem is ill-posed. However, one anticipates that the stability will increase 
if one increases the frequency. This phenomenon was observed numerically in the 
inverse obstacle scattering problem |l5l. Several rigorous justifications of the in- 
creasing stability phenomena in different settings were obtained by Isakov et al 
[l6l|7l[8l[T0l[IIl. Especially, in |[8l, Isakov considered the Helmholtz equation with 
a potential 

- An - A;^n + gn = in VL. (1.1) 
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He obtained stability estimates of determining q by the Dirichlet-to-Neumann map 
for different ranges of A;'s. All of these results demonstrate the increasing stability 
phenomena in k. For the case of the inverse source problem for Helmholtz equa- 
tion and an homogeneous background it was shown in Q that the ill-posedness 
of the inverse problem decreases as the frequency increases. 

In this paper, we study the acoustic wave equation. Let Vt C W be a bounded 
domain, where n > 3. Let dVL be smooth. We consider the equation 

(A + A;2g(x))u(x) =0 in Vt, (1.2) 

where the real-valued q{x) is the refractive index. Assume that the kernel of the 
operator A + k'^q(x) on Hq{^1) is trivial. Associated with (11.21) . we define the 
Dirichlet-to-Neumann map (DN map) A : H^/'^{dVt) H-^/'^{dVt) by 



do. 



where u is the solution to (11.21) with the Dirichlet condition n = / on dVt, and u 
is the unit outer normal vector of dVt. The uniqueness of this inverse problem is 
well known llT3l . This inverse problem is notoriously ill-posed. For this aspect, 
Alessandrini proved that the stability estimate for this problem is of log type flTl 
and Mandache showed that the log type stability is optimal In this paper, we 
would like to focus on how the stability behaves when the frequency k increases. 
Now we state the main result. 

Theorem 1.1. Assume that qi{x) and q2{x) are two sound speeds with associ- 
ated DN maps Ai and A2, respectively. Let s > (n/2) + 1, M > 0. Suppose 
\\(Ii\\h^{q.) < M {I = 1,2) and supp(gi — ^2) C Q. Denote q a zero extension of 
Qi — Q2- Then there exists a constant Ci, depending only on n, s, and Q, such that 
ifk"^ > l/(CiM) and ||Ai - A2II* < 1/e then 



C ( \ 
||glliT-(R") < y^exp(CA:2)||Ai - A2II, + C P + log — 

(1.3) 

holds, where (7 > depends only on n, s, ^2, M and supp(gi — q^). Here ||-||* is 
the operator norm from H^^'^{dil) into H~^^'^{dVt). 

Remark 1.2. 1 . The estimate (11.31 ) consists two parts - Lipschitz and logarithmic 
estimates. As k increases, the logarithmic part decreases and the Lipschitz part 



-{2s-n) 
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becomes dominated. In other words, the ill-posedness is alleviated when k is 
large. 

2. We would like to remark on the constant C exp{Ck'^) / k'^ appearing in the Lips- 
chitz part of (11.31) . 1/k'^ comes from k'^q in the equation, which appears naturally, 
while, exp(CA;^) is due to the fact that we use the complex geometrical optics so- 
lutions in the proof. Even so, we expect that the there is an exponential growth of 
the constant with frequency since we do not assume any geometrical restriction 
on q{x) other than regularity. For the wave equation it has been shown by Burq 
for the obstacle problem p| that the local energy decay is log-slow and this is due 
to the presence of trapped rays. Notice that in our case we can have trapped rays. 
For the case of simple sound speeds we expect that there is no exponential increase 
in the constant. In W2\ a Holder stability estimate was obtained for the hyperbolic 
DN map for generic simple metrics. For very general metrics there is not known 
modulus of continuity for the hyperbolic DN map, see for convergence results. 

However, in practice, k is fixed and so is the constant. Therefore, one should 
expect to obtain a better resolution of q from boundary measurements when the 
chosen k is large. 

3. Unlike the result in [[8l Theorem 2.1] (for equation (11.11 )) where the stability 
estimates were derived in different ranges ofk, estimate (11.31) is valid for all range 
ofk provided k^ > l/((7iM) . 

The proof of Theorem 11.11 makes use of Alessandrini's arguments [HI and the 
CGO solutions constructed in [[T3l . The main task is to keep track of how k 
appears in the proof of the stability estimates. 

2 Complex geometrical optics solutions 

In this section, we construct CGO solutions to the equation (11.21 ) by using the idea 
in |[T3l . The main point is to express the dependence of constants on k explicitly. 
We first state two easy consequences from the results in ffT3l . 

Lemma 2.1 (see |fT3l Proposition 2.1 and Corollary 2.2]). Let s >0 be an integer 
Let £q > 0. Let ^ G satisfy ^ ■ ^ = and \^\ > Eq. Then for any f G H'^IQ) 
there exists w G H'^{Q) such that w is a solution to 

Aw + ^ - Vw = f inQ 
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and satisfies the estimate 



\w\\h-{q,) < 1^11/ linden), 



where a positive constant Cq depends only on n, s, Eq and Q. 
By using this lemma, we can obtain a solution to the equation 

A^ + e-V^A + # = / (2.1) 
satisfying some decaying property as in the following lemma. 

Lemma 2.2 ( lfT3l Theorem 2.3 and Corollary 2.4]). Let s > n/2 be an integer. 

Let Eq > 0. Let ^ G C" satisfy ^ ■ ^ = and \^\ > Eq. Let f,g e H'{Q). Then 
there exists Ci > depending only on n, s, Eo and Q such that if 

\^\>C^\\g\\Hsin) 

then there exists a solution ip G H'^{Q) to the equation (12.11 ) satisfying the estimate 

where Cq is the positive constant in Lemma \2.1\ 

The needed CGO solutions are constructed as follows. 

Proposition 2.3. Let s > n/2 be an integer Let Eq > 0. Let ^ E C" satisfy 
,^ ■ ,^ = and \^\ > Eq. Define the constants Co and Ci as in Lemma IZ2l Then if 

\^\>Cie\\q\\H^(n) 
then there exists a solution u to the equation (11.21) with the form of 



u[x] 



exp(^^-x^ (l + ^(x)), (2.2) 



where ijj has the estimate 



H={n) S — — \\(1\\h''{q)- 



Proof Substituting (|2.2I) into (11.21) . we have 

A^^ + ^ ■ + k^qiJ = -k^q. 
Then by Lemma [T2l we obtain this proposition. □ 
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3 Proof of stability estimate 

This section is devoted to the proof of Theorem ll.il We begin with Alessandrini's 
identity. 

Proposition 3.1. Let ui be a solution to (|1.2I) with q = qi, then we have 

k'^ / {(l2- qi)uiU2dx = {{Ai- K2)ui\an, U2\avt)- 

Now we would like to estimate the Fourier transform of the difference of two 
g's. We denote J-'{f ) the Fourier transformation of a function /. 

Lemma 3.2. Let s > {n/2) + 1 be an integer and M > 0. Assume WqiWu^CD.) < 
M, supp(gi — ^2) C and > 1/CiM, where Ci is the constant defined in 
Lemma |272] corresponding to Eq = 1. Let q be a zero extension of qi — q2 and 
do > Ci. Suppose that X £ C^{fl) satisfies x = 1 f^^^^ supp(gi — ^2)- Then for 
r > and rj G M" with \ri\ = 1 the following statements hold: ifO < r < aok'^M 
then 

i II II i 

\Tq{rri)\ < "^'""^^"^ ||g||H-.(M.) + - exp(CaoA;2M)||Ai - A2II* (3.1) 
holds; ifr> Cik'^M then 

\:Fq{rT^)\ < + gexp(Cr)||Ai - A2II. (3.2) 

T K 

holds, where C > depends only on n, s and fl. 

Proof. In the following proof, the letter C stands for a general constant depending 
only on n, s and O. By Proposition [23l we can construct CGO solutions ui{x) to 
the equation (11.21) with q = qi having the form of 

ui{x) = exp ■ (1 + ipiix)) 

for / = 1, 2, and we have 

(g2 - qi) exp Q(6 + 6) ■ a;^ (1 + ^1 + ?/;2 + i^ii'2) dx 

- ^2)ul\aQ., U2\dn) (3.3) 
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from Proposition [3T1 where ipi satisfies 

Ck'^ 

\\'ipi\\H''{n) < -rpr\\qi\\Hs{n) 

if^i e C" satisfies ■ = 0, > I and 

161 > C,e\\qi\\Hsin). (3.4) 
We remark that < C also holds. Indeed, we have 

„ „ Ck^\\q,\\ 

"*"«•'"' - — \U\ — - c,p||„i|h.,«, ^c'r'^- 

Now, let r > 0, and ?7 G R" satisfy |?7| = 1. We assume that a, C G M" satisfy 
a-r] = a-C = r]-C = Oand |Cp = |ap + r^. (3.5) 
Define .^i and ,^2 as 

= ( + ia — irr] and .^2 = — C — i<y — irrj. 

Then we have 

6 ■ 6 = 0, = ICr + |«P + = 2|CP (/ = 1, 2) and 1(^1 + 6) = -^rr/. 
Hence by (13.31 ). we immediately obtain that 

J^q{rr]) = - / (^2 - ^i) exp(-ir?7 ■ x){iJi + + i'ii'2) dx 
Jn 

+ ^(^^^ ~ A2)ni|aQ, ^219^) (3.6) 
provided |.^/| > 1 and (|3.4I) are satisfied. We first estimate the first term on the 
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right hand side of (13.61) by 



(^2 - qi) exp{-irr] ■ a;)(V'i +4^2 + V'iV'2) dx 



(92 - qi) exp(-zr?7 ■ + ^2 + ^1^2) c^a^ 



< 11^2 - + V^2 + V^iV^2)||^.(^) 

< ||^|H-'=(R")llx||//=(f^)(il^ilk-{n) + ||^2||H=(f7) + \\i^i\\H'^in)\\i^2\\H'^[n)) 



< llfll«-.,..,llx-ll«.,«, + ^ + c:^) Ell«ll« 



(a) 



ICI 



||gllif-=(M")Xlll*ll^^' 



1=1 



since + "^2 + '^i'^2) ^ -^o(^) ^'^d > ''^Z^- 

On the other hand, by taking R large enough such that f2 C -Bij(O), we have 



< C exp 

Likewise, we can get that 
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i? (1 + C) = Cexp 



||Vn,| 



an 



^u, + exp 



(VV'/) 



< CICIexp (^Mi?^ +Cexp (^^i?^ HVMns-^^n) 

< CiClexp (^Mi?^ +Cexp (^^i?^ UnllHsin) 
<Cexp(C|C|) 
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since s — 1 > n/2. Consequently, we have 

<Cexp(C|C|). 

Therefore, we can estimate the second term of the right-hand side of (13.61) by 

(A1-A2) 

<Cexp(C|C|)||Ai-A2|U. 

Summing up, we have shown that for r > and for 77 G M" with |?7| = 1 if we 
take a and C, satisfying the conditions (|3.5I) . |C| > 2^^/^ and 

\C\>2~^'^Cie\\qi\\Hsin) (3.7) 

then 

+ ^exp(C|C|)||Ai-A2|U (3.8) 

holds. 

Now assume that < M and k"^ > 1/CiM. Thus if 

Id > Cik^M (3.9) 

holds, then (13.71 ) and \(\ > 2~^/^ are satisfied. Pick oq > Ci. We first consider 
the case where < r < a^k'^M. By choosing a and ( satisfying 

a-n = a-C = ri-C = 0, \C\ = aok'^M{> r) and \a\ = ^/{aok^My - 

both (13.51) and (13.91) are then satisfied since ao > Ci. Hence we obtain (|3.8I) . that 
is (13.11) . On the other hand, when r > Cyk'^M, we can choose a = 0, 77 ■ C = 
and Id = r. Then (|3.5I) . (13.91 ) are satisfied and thus (13.81) holds and consequently 
(1X21) is valid. □ 

Now we prove our main result. 
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Proof. As above, C denotes a general constant depending only on n, s and Vt. 
Written in polar coordinates, we have 



= c 



\J^q{rr])\\l + r')-'r'^-' dr]dr 



JH=i 

aok^M 



+ [ [ \J^q{rr])Wl + r^y'r''"Ur]dr 

Jaok'^M J\ri\=l 

+ f [ \Tq{r7])Wl+rY'r''-Ur]dr] 
Jt J\7i\=1 J 



\J'q{rri)\'^{l + r'') V" ^drjdr 



r,\=l 



--: Cih + h + h), 



(3.10) 



where ao > Ci and T > aok'^M are parameters which will be chosen later. Here 
Ci is the constant given in Lemma [X2l From now on, we take k"^ >!/ (CiM). 
Our task now is to estimate each integral separately. We begin with I^. Since 

\J^q{rr])\ < C\\qi - q2\\L^n), Qi - q2 ^ H^i^), and s > n/2, we have that 



T 



ki - 92||i.(n)(l + r^rr--' dr < CT--\\q, - q^Wh^^^ 



< CT " ( e\\qi - q2\\H-sin) + " mH^{n) 



-I 

e 



< CT-^ e\\q\\'^. + — 



for e > 0, where m := 2s — n. 

On the other hand, by Lemma [3]2l we can estimate 



(3.11) 



aok^M 



'1 +r')-V"-Mr 



X 



2 ll^l/i'-''(R") 



, exp(2Caofc^M) „^ 
H ~ U^-i ~ ^1-2 



< C I (1 + r^)-'r''-Ur 



a 



X ii;ri|2 



ll^l 



Cexp(CaoPM) 



Ai-A2||^, 



(3.12) 
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where x ^ C'cr(^) satisfies x = 1 near supp(g2 — 9i) and := \\x\\H^[n)- In 
view of 

(1 + r^)-^""^ dr<[ dr < C {a^k'^ M)-^'+'^-'^ 

and 

-T 



POO 

<exp((7r) / (1 + r2)-V"-idr 
Jo 

< Cexp(CT), 

we have that 

Jaok'^M 

+ ^||Ai-A2||^ r exp(Cr)(l + -2^--"-i 



r')-'r''-'dr 

aok^M 



cc"^ c 

< -^tMh-h^-) + Tiexp(Cr)||Ai - A^ll^ (3.13) 



Combining (IXTOl )- (l3J3]) gives 

||gllH-.(M.) < C(/i + /2 + /3: 

2 



U;Q A/ 

+ -/II^Ih-^(M") + Tiexp(CT)||Ai - A^ll^ 

"o / 

C CM"^ 
+ ^(exp(Caofc'M) + exp(Cr)) ||Ai - K2\\l + 
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where positive constants C2 and C3 depend only on n, s and Vt. 
Now we pick and e as 

ao = - C*! and e = — - 

(if needed, we take C2 large enough). We then obtain that 



^exp{Cak'^)A + C^{T) (3.14) 

rC 



for T > aoA;2M = 2C2C^k'^M = ak'^, where 

$(T) := exp(C4r)A + M^T-^"^, 

rv 

A := ||Ai - A2||^, a := 2C2C^M'^ and C4 > depends only on n, s and 
To continue, we consider two cases: 

ak^<p\ogj (3.15) 

and 

aP>p\ogj, (3.16) 

where p will be determined later (see (13.241) ). 

For the first case (|3.15l) . our aim is to show that there exists T > ak'^ such that 

$(T)<2C5(^A;2 + logi^ . (3.17) 

Substituting (13.171 ) into (13.141) clearly implies (fT31) . Now to derive (13.171 ), it is 
enough to prove that 

1 / . \ -2m 

-exp(C4T)A < C5 + log-j (3.18) 
and ^ 

11 



Remark that (13.191) in equivalent to 
which holds if 

T > Cg-^/^^M^/™ (l + ^) log ^ (3.20) 

because of (I3l31) . Setting T = plog(l/A) (> ak"^ by (irBl) ). then (Il20l) holds 
provided 

Now we turn to (|3.18|) . It is clear that (|3.18l) is equivalent to 

Cap log ^ < log C5 + 2 log /c^ + log ^ - 2m log [k"^ + log ^ j (3.22) 
since T = plog(l/A). It follows from (13.151 ) that 

log ^/c^ + log < log log ^ + log ^ = log + + log log ^. 
Hence (13.221 ) is verified if we can show that 

CiP log ^ < log C5 - 2 log(MCi) + log ^ - 2m (log f ^ + l) + log log \ 
A A \ \a / A 



I.e. 



(1 - C4P) log ^ - 2m log log ^ + log C5 - 2 log(MCi) - 2m log (^^ + l) > 

(3.23) 

for log(l/A) > 1. Now we choose 



Then (13.231 ) becomes 

log4-4mloglog4 + 21ogC5-41og(MCi)-4mlogf- + l) > 0. (3.25) 
A A \a / 
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Notice that 



inf 1 loff — — 4m loff log — I = inffz — 4mlo£rz) 

0<A<l/e \ ^ A A) z>l^ ^ ' 



g 

> inf (2 — 4m log 2;) = 4m log . 

2>o' ' 4m 

Hence if we choose C5 such that 

C5>(MCi)^(^ + l)""(^y'" (3.26) 
then (13.251) follows. Finally, we take 

C5 := max \c\ , p'^-j (l + ^)'" , 

which depends only on n, fi, s, M and x- With such choice of C5, the conditions 
(13.261 ) and (13.211) hold, and thus estimate (13.171 ) is satisfied. 

Next we consider the second case (13.161 ). By (13.141 ) with T = a/c^, we get that 

||g|lH-»(Rn) < ^ exp{Cak'')A + ^ exp{C^ae)A + CM2(aA:2)-2- 



Hence it remains to show that 

fc-^™ < a ( A;2 + loe 



\ —2m 



I.e. 

e>c,'"'-[k' + \og^. 021) 

Since 

fc^ + logl< (1 + ^) 
by (13.161) . we have (13.271 ) if we take large enough so that 

2m 

Ce > 1 + - 
P. 

The proof is completed. □ 
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